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We report strongly momentum-dependent local charge screening dynamics in CE-type charge, 
orbital, and spin ordered Lao.sSri.sMnO,!, based on Mn K-edge resonant inelastic x-ray scattering 
data. Through a comparison with theoretical calculations, we show that the observed momentum 
dependence reflects highly localized, nearest neighbor screening of the transient local charge per- 
turbation in this compound with an exciton-like screening cloud, rather than delocalized screening. 
The size of the screening cloud is estmated to be about 0.4-0.5 interatomic distances. 

PACS numbers: 78.70.Ck, 78.70.Dm, 71.27.+a, 75.47.Lx 



The dynamic screening of the Coulomb interaction 
plays a central role in determining the electronic proper- 
ties of materials pQ . The response of valence electrons to 
a potential, in particular on time scales of the order of 
femtoseconds, is through excitation of electron-hole pairs 
which screen "bare" charges in the system. The screen- 
ing is described theoretically by the density-density cor- 
relation function, or its Fourier transform, the dynamic 
structure factor [2]. Spectroscopies that probe the va- 
lence band, such as photoemission, are sensitive to these 
screening dynamics. However, it is often not obvious how 
to separate the kinetics of a charged particle and the re- 
sponse of the rest of the system to its presence [3]. In 
contrast, core- level spectroscopies provide an alternative 
way of studying the screening dynamics. By removing 
a deep-lying core electron, a strong local potential is 
created that exists for a very short time, i.e., the core 
hole lifetime. Essentially, one creates a short-lived lo- 
calized "test" particle, and measures the response of the 
electrons to this local transient potential. This type of 
screening dynamics has drawn the attention of scientists 
for decades @HS] ■ 

For transition metal compounds, K-edge resonant in- 
elastic x-ray scattering (RIXS) [7 offers the intriguing 
possibility of projecting the excitations related to the core 
hole screening onto valence band excitations. Specifically, 
it has been shown that K-edge RIXS can be directly re- 
lated to the dynamic structure factor in the limit of a 
strong or weak core hole potential, U corc , relative to the 
band width [8] . For the case where U core is comparable to 
the bandwidth, more typical for 3c? transition metal com- 
pounds, the screening is more complicated because there 



is an asymmetry between the electron and hole excita- 
tions and the intermediate states can not be integrated 
out In this case, the RIXS response is believed to be 
sensitive to the transient screening of the intermediate 
states to the core hole potential [5]. 

In this letter, we report RIXS measurements of the mo- 
mentum and energy dependence of the screening dynam- 
ics for a transient local potential in a CE-type charge, 
orbital, and spin ordered manganite, Lao.sSri.sMnO^ 
We find strong momentum dependence of the intensity 
of the across-gap excitation, with a dramatic increase on 
moving away from the two dimensional (2D) zone center. 
We show that this behavior reflects the size and shape 
of the real-space screening cloud and demonstrate that 
in Lao.5Sri.5Mn04, the screening distance is very short, 
with a screening cloud of about 0.4-0.5 interatomic dis- 
tances in size. 

A single crystal of Lao.5Sri.5Mn04 was grown by the 
traveling solvent floating zone method. It has a tetrag- 
onal structure at room temperature with H/mmm sym- 
metry and undergoes a charge and orbital ordering tran- 
sition around 230 K, accompanied by complex structural 
distortions [10) . For simplicity, we use here the H/mmm 
notation throughout. The wave vectors of the charge and 
orbital ordering are then of the form (|, |, L) and (|, \, 
L), respectively. In the low temperature ordered state, 
Lao.5Sri.5Mn04 is an insulator with a large gap between 
the e g states [TT]. The behavior of the excitation be- 
tween these predominantly Mn 3c? states, labeled as a 
d-d transition, is the focus of this study. The Mn K-edge 
RIXS experiments were performed at Advanced Photon 
Source on beamlines 30-ID and 9-ID with an instrumen- 
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tal energy resolution of about 270 meV (FWHM). The 
polarization dependence of the RIXS process is controlled 
by placing the [001] and [110] directions of the crystal in 
the scattering plane. The incident beam polarization is 
perpendicular to the scattering plane, i.e., parallel to the 
[110] direction. Thus the incident polarization condition 
is fixed for all the Q = (H, H, L) points surveyed. All the 
data presented were collected at T = 20 K, well below 
the Neel temperature(llOK) [T^]. Data are normalized 
by incident beam intensity and corrected for footprint 
variations. 

In Fig.JIJa), we show RIXS spectra taken at three Q 
points. In each case, there is a large elastic line cen- 
tered at zero energy loss. The d-d transition appears 
as a peak on the tail of the elastic scattering at around 
2 eV, consistent with optics [T3] and EELS [H] obser- 
vations. This across-gap transition has also been ob- 
served by K-cdgc RIXS on other manganites [HI ITo] . 
Remarkably, the RIXS spectra show a strong momen- 
tum dependence of the intensity of this feature. At 
Qo = (—0.03, —0.03, 7.20) with very small in-plane mo- 
mentum transfer, the 2 eV peak is almost unobservable. 
This momentum dependence is confirmed with RIXS 
spectra collected at a large number of Q points, as shown 
in Fig.Qb) andjljc) . To control the systematics resulting 
from polarization factors, the data were taken either at 
fixed sample angle, 8, or fixed detector angle, 28. These 
conditions result in data taken along three lines in re- 
ciprocal space. In all cases, the incident polarization is 
parallel to the [110] direction. With the detector position 
28 fixed, polarization effects associated with the outgo- 
ing x-ray are eliminated. Figure[ljc) shows RIXS spectra 
with the elastic line subtracted [T7] for the Q points with 
28 = 68°. We take the integrated intensity, I(Q), over 
the 1-3 eV range as a measure of the strength of the 2 eV 
peak. The size of the symbols in Fig.jljb) is proportional 
to /(Q). A clear systematic dependence on momentum 
transfer is observed. 

The integrated intensities of the d-d excitation are plot- 
ted as a function of the in-plane momentum transfer in 
Fig.[2jb). In order to quantitatively compare the ex- 
perimental data with the theoretical calculations (dis- 
cussed below), the integrated intensities are plotted rel- 
ative to the intensity at Q = (-0.03,-0.03,7.20) [15] . 
i.e., /(Q) — /(Qo). This removes the uncertainty in de- 
termining the common background for all Q points. The 
strength of the 2 eV d-d excitation exhibits a minimum 
at zero in-plane momentum transfer and a maximum at 
(0.5, 0.5, L). Interestingly, although there is a large 
variation in the L values for the various Q points [see 
Fig.jljb)], all the measurements collapse onto a single 
curve in Fig.[2](b). This demonstrates that there is neg- 
ligible L dependence to this behavior, a result consis- 
tent with the 2D nature of this single layered manganite. 
Further, it implies that the polarization factors are in- 
deed constant for our experimental geometry. From here 
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FIG. 1. (a) RIXS spectra at three Q points, (b) The Q 
points surveyed in the (H, H, L) plane. The radius of the 
dot is proportional to the integrated intensity of the 2 eV 
peaks. 9 and 29 are the incident and detector angles, (c) 
RIXS spectra for the Q points along the 29 = 68° line, with 
the elastic intensity subtracted. The grey shaded region is the 
energy window used in calculating the integrated intensity of 
the feature. 



on, the momentum transfer will be denoted simply as 
Q2D = {H,H) since the L component is irrelevant. 

The experimental data in Fig.[T] and Fig. [2] show our 
main experimental observations. The across-gap d-d ex- 
citation in Lao.5Sri. 5 Mn04, as observed in the RIXS pro- 
cess, exhibits a strong momentum dependence. While 
the position of the peak shows no appreciable dispersion, 
the intensity increases rapidly as the in-plane momentum 
transfer increases away from the 2D zone center. Near 
the zone center, the spectral weight of the 2 eV feature 
almost disappears. This is a surprising result. The dis- 
appearance of this RIXS spectral weight at the 2D zone 
center cannot be the result of the dynamic structure fac- 
tor going to zero, since this feature is still observed in 
the optical response |13j . which probes the zero momen- 
tum transfer response function. This demonstrates that 
K-edge RIXS in Larj.sSri.sMnC^ is indeed in the inter- 
mediate core-hole potential regime, discussed in the in- 
troduction. In the following, we detail this momentum- 
dependence and show that it arises from the intermediate 
state screening dynamics and in particular that it reflects 
the real-space extent of the screening cloud. 
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To understand this strong in-plane momentum de- 
pendence, we calculated the RIXS response from 
Lao.sSri.5Mn04 for a two-dimensional 16 x 16 Mn cluster 
with periodic boundary conditions. The initial and final 
states of the unperturbed system, and the intermediate 
states in the presence of the Is core hole on-site Coulomb 
potential, are solved numerically with a tight-binding ap- 
proach. The Hamiltonian employed is similar to the one 
in Ref. [T5], which includes the nearest-neighbor elec- 
tron hopping within the Mn02 plane, the Jahn- Teller 
and isotropic electron-lattice coupling, the Hund's cou- 
pling to the CE-type ordered t 2g spins, and the Coulomb 
interaction between e g electrons within the Hartree-Fock 
approximation 20 . The sizes of the distortions of the 
oxygen octahedra are taken from Ref. [TU]. The RIXS 
spectra are then calculated from the Kramers-Heisenberg 
formula [7112]: 



/ 



E 



(f\V*\n)(n\V\g) 



E„ + Huju - E n + iTr. 



8{E f -E g - HAuj), 
(1) 

where |/), |n), and \g) represent the final, intermediate 
and initial states, and Ej, E n and E g their energies. T n 
is the inverse of the intermediate state lifetime, and T>'^ 
and T> are the RIXS dipole transition operators. Tiwk and 
hAiv are the incident x-ray energy and the energy loss, 
respectively. The calculated RIXS intensity is averaged 
over configurations in which the zig-zag chains of orbital 
order are along either the [110] or the [110] directions, to 
take into account twining effects in real crystals. Details 
of the calculation will be published elsewhere. 

The calculated RIXS spectra were found to be most 
sensitive to the e g -e g hybridization and the coupling of 
the e g electrons to the distortions of oxygen octahedra. 
These two effects are parametrized as to and A in our 
Hamiltonian, where to is the hopping between 3x 2 -r 2 
orbitals along the x direction and A is proportional to 
the strength of the clectron-phonon coupling [20 . With 
reasonable parameter values 20 and the combination of 
to = 0.9 eV and A = 7.41 eV/A, the calculated spectra 
shown in Fig.[2]ja) and the thick (blue) line in Fig.[2]jb) 
closely resemble the experimental observations. The in- 
tensity of the calculated RIXS response peaks near 2 eV, 
and increases rapidly as Q 2 d increases away from (0, 0), 
towards (0.5, 0.5), as seen in the experiments. We note 
that the calculated spectra in Fig.[2ja) suggest a slight 
dispersion of about 130 meV of the 2 eV peak, which is 
much smaller than that reported for LaSi^Mj^Oy [IB] , 
Such a small dispersion, roughly equal to the experimen- 
tal step size taken in Fig.JlJc), is below the detection limit 
of our experiment. 

The sensitivity of the RIXS response to intersite 
hopping and the electron-phonon coupling is shown in 
Fig.[2jb) by varying to and A [2U]. For a given t , A is 
constrained such that the d-d excitation in the RIXS re- 




FIG. 2. (Color online) (a) Contour plot of RIXS intensity 
calculated for the electron hopping parameter to = 0.9 eV, 
after averaging over twin domains, (b) The integrated RIXS 
intensity of the 2 eV peak relative to the (-0.03,-0.03) point, 
plotted with respect to the in-plane momentum transfer along 
the (H, H) direction. Symbols represent experimental data. 
Lines represent theoretical results for different values of to. 
Both experimental data and theoretical results are normalized 
for comparison. 



sponse peaks near 2 eV. Henceforth, only to is mentioned 
for simplicity. The details of the combinations of to and 
A can be found in the supplemental material [20] . As was 
done for the experimental data, the calculated response 
is integrated over the same 1-3 eV window to generate 
the curves in Fig.[2jb), and again the value at Qo =(- 
0.03, -0.03) is subtracted. The calculated results show 
the best agreement with the experimental observations 
when to = 0.9 eV. For larger to values, the calculated 
RIXS response differs significantly from the experimen- 
tal data. Thus our study sets the upper limit of to- The 
inability to precisely determine the parameter values is 
largely due to the difficulty in determining the contribu- 
tions from other inelastic scattering processes that give 
rise to a smooth "background" in the low energy loss 
region. 

To further understand the implications of the observed 
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FIG. 3. (Color online) (a) and (b): Screening configuration 
in real space for to = 0.9 and 1.5 eV, respectively. The top 
left half corresponds to the case with the core hole at a Mn 3+ 
site, while the bottom right it is at a Mn 4+ site. The volumes 
of the red and blue spheres are proportional to the electron 
and hole numbers. The big red spheres at Mn 4+ core hole 
sites represent about 0.9 electrons, (c) and (d): Integrated 
RIXS intensity plotted in the (H, K, 0) plane of reciprocal 
space for to = 0.9 eV and to = 1.5 eV, respectively. Red 
and blue represent the maximum and minimum intensities, 
respectively. 



momentum dependence of the RIXS spectrum for the 
screening dynamics, we calculate the real-space screen- 
ing configurations from the lowest energy intermediate 
eigenstates. These are shown in the top panels of Fig. [3] 
Figures [3^a) and[3jb) compare the charge redistributions 
for to = 0.9 and 1.5 eV, respectively, with the core hole 
at either a Mn 3+ or a Mn 4+ site. The volumes of the red 
and blue spheres scale with the screening electron and 
hole densities on the individual sites. Figures |3jc) and 
J3jd) show the calculated RIXS intensities for the two t 
values over half of a Brillouin zone. 

As expected, the excited hole distributions are more 
localized near the core hole sites for the smaller value of 
to. For to = 0.9 eV, the screening hole is tightly bound 
to the excited electron with more than 90 % of the ex- 
cited charge located on the three nearest neighbor sites 
along the zig-zag chain. The predominant wave vector for 
these electron and hole distributions is (0.5, 0.5), coinci- 
dent with the location of the RIXS peak intensity max- 
imum in reciprocal space in Fig.[3]jc) and in agreement 
with our experimental results. For the large hopping pa- 
rameter, to — 1.5 eV, the screening pattern in real space 
changes drastically. The majority of the hole distribution 
in Fig.[3][b) is now beyond the nearest neighbor sites, and 
is spread throughout the zig-zag chains. This difference 
in screening dynamics is directly reflected in the RIXS 
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FIG. 4. (Color online) The excited hole number normalized 
by the excited electron number at the core hole site, plotted 
in semi-logarithmic scale with respect to the distance from 
the core hole site along the zig-zag chain. 



response, with the maximum of the RIXS response then 
shifted to around (0.25, 0.25), as shown in Fig. [2^b) and 
Fig.|3jd). This pattern is completely at odds with that 
seen in the experiment results. 

The relationship between the hopping strength and the 
charge redistribution in the screening process apparent 
in the theoretical calculations is shown more clearly in 
Fig.|4j where the relative hole number at a given site 
is plotted as a function of the distance from the core 
hole site, as measured along the zig-zag chain. The semi- 
logarithmic plot reveals an exponential decay of the hole 
density, confirming the presence of exciton-like screening 
clouds. The size of the screening cloud, which charac- 
terizes the screening dynamics and determines the RIXS 
response, depends strongly on the hopping strength to- 
Taking the to — 0.9 eV case, which best describes the 
RIXS data, we fit the excited hole distributions to an ex- 
ponential function and find that the size of the screening 
clouds are 0.4 and 0.5 atomic spacings for the Mn 3+ and 
Mn 4+ sites, respectively. 

In summary, we observe a highly momentum- 
dependent K-edge resonant inelastic x-ray scattering 
intensity in the orbital ordered, layered manganite 
Lao.sSri.sMnC^. We interpret this through a comparison 
with calculations based on a tight-binding approach, and 
show that these observations imply a highly localized, 
nearest neighbor screening of the local charge perturba- 
tion. We further find that the momentum dependence 
of the RIXS spectrum reflects the pattern and range of 
the screening in real space, and thus are able to measure 
the size and shape of the screening cloud. We find that 
the screening cloud is localized to a few Mn sites in the 
Mn-0 plane, emphasizing the short range nature of the 
Coulomb interactions in the manganites. 
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